On the Ginzburg-Landau Free Energy of Charge Density Waves 
with a Three-Dimensional Order 
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The effective free energy of a charge density wave (CDW) with a three-dimensional order is derived 
from a microscopic model (Frolich model) based on the path integral method. Electron hoping and 
Coulomb interaction between chains are taken into account perturbatively, leading to an elastic 
interchain coupling of the CDW ordered state. 



I. INTRODUCTION 

After a pioneering study by Peierls , charge density 
wave (CDW) has been a subject of continuous interest 
because of its unique nature as an electronic condensate. 
Especially, its behavior under an electric field was inves- 
tigated by Frolich as a possible origin of superconduc- 
tivity. Although it turned out that the CDW mechanism 
does not apply to the superconductivity, many interest- 
ing behaviors have been found in the electromagnetic 
response of CDW Even today, experimental de- 

velopments are adding intriguing new phenomena to the 
book of CDW. For example, AB effect in CDW's with 
columnar defects || , field effect [|| and current effect |t]] 
of CDW transport have been discovered although some 
theoretical developments are needed for the full under- 
standing of these phenomena. 

In a preceding paper [|| , we have studied CDW's in a 
transverse (perpendicular to the chains) electric field and 
pointed out a possibility of a state similar to the mixed 
state of superconductors. Although interchain coupling 
plays an essential role in this theory, it was only phe- 
nomenologically introduced. The aim of the present pa- 
per is to give a microscopic foundation of the interchain 
coupling by deriving it from the electron-phonon model. 

Although there have been proposed several origins 
of the interchain coupling |s|Jic(|, the derivation of the 
Ginzburg-Landau (GL) type effective free energy taking 
them into account in a unified framework has not been 
given until now. We reexamine the effective free energy 
of the CDW taking account of the electron hopping and 
the Coulomb interaction between chains as sources of the 
interchain coupling. These effects are treated pertur- 
batively employing the path integral formulation. It is 
shown that the free energy previously used by the present 
authors || can be derived from a microscopic model. 

For simplicity we limit our discussion to the case with- 
out impurity and commensurability pinning. 



II. FREE ENERGY OF A CDW SYSTEM 

We consider a three-dimensional (3-D) stack of f-D 
conducting chains, each of which is described by Frolich 



model (coupled electron-phonon system) . At low temper- 
atures, this system is expected to show a 3-D CDW order 
due to interchain coupling P,§|. We treat the Coulomb 
interaction and the electron hopping between chains as 
origins of interchain coupling. However the electron hop- 
ping is assumed to be small so that it does not break the 
one-dimensionality of the system |pT|| . 



by. 



A. One-dimensional electron-phonon system 



The action of the electron-phonon system S e . p is given 



Sc-p — 



Jo 



9 



(2.1) 



where i is the index numbering the chains, L is the length 
of the chains and, m and — e are the mass and charge of an 
electron, respectively. The chemical potential fi is given 
by fi = h 2 kp/(2m) with kp being the Fermi wave num- 
ber. We assume that /i is common to all the chains |lT| . 
The phonon frequency and electron-phonon coupling con- 
stant of wave number q are denoted by uj q and g q , re- 
spectively. Electron field on the i-th chain with spin a is 
represented by ipi a (x,T) and its Fourier transformation 
is given by a,i a {k,T) = L -1 / 2 j dxtpi a (x,T)exp(—ikx). 
Note that the summation over a is suppressed in the 
present paper. The Fourier component of phonon field of 
the i-th chain is denoted by bi(q) which is related to the 
displacement of ions Ui(x) by, 



Ui{x) = -= ^2 a q e q {btiq) + b*(-q)} 



(2.2) 



1 



where a q = %/ \j2pMfwjq with p M being the average 
mass density of ions and e q — qj\q\. Here we have consid- 
ered only the phonons corresponding to the displacement 
of ions parallel to the chains. 

We consider the Peierls instability at the wave number 
Q = ±2kp and treat only the phonons with wave number 
close to ±2fcp. It is convenient to divide ipi a into two 
parts as ipi a = e lkFX Ri a + e~ tkFX L ia , where Ri a and L icr 
stand for right and left moving electrons, respectively, 
and are assumed to be slowly varying fields. 

The complex order parameter (or energy gap) of CDW 
is introduced by 



Ai(x) = -j= 5>*(-Q -q) + bi(Q + q))i 



iqx 



(2.3) 



where (• • •) denotes the statistical average and we neglect 
the g-dependence of gQ+ q , writing it as g. Here we have 
assumed that the phonon field and the order parameter 
are classical fields which do not depend on imaginary 
time r. Then the action can be rewritten as, 



"Sc-p — 



where 
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dx 



2g 2 



fc< +) A,(x) 
(-) 



A* (or) k, 
L ia (x,T) 



(2.5) 
p. (2.6) 

Here g-dependence of 0Jq+ q is ignored as ujQ+ q — > ujq. 
By applying the variation with respect to A*(x) we ob- 
tain self-consistent equation, 



hd T (±ihp + d x Y 



fiu>Q 



(L*a(x,T)R ta {x,T)}. 



(2.7) 



where (fi(x, t) is the value of the scalar potential ip(r, t) 
on the «-th chain and pi(x,r) is the total charge den- 
sity on the i-th chain. We write pi (x,t) = pf~(x,r) + 
p 1 ° n (x, r), where pf(x, t) and p 1 ° n (x, r) are the electronic 
and the ionic contribution to the charge density, respec- 
tively. 

The electronic contribution pf has two characteristic 
parts, the long range (q ~ 0) and the short range (q ~ 
2kp) charge modulation, which we denote by Pq\ and p < Q i , 
respectively. These quantities are written as, 



Pi = Po 



POt — ~ e {R*<jRi< 

4, = e {LlR^e*** + R* a L i(7 e~ iQx ) 



(2.9) 



The ionic contribution pf n also has two contributions 
Pg™ and Pq^. For ions Pq°" is a constant, since we take 
account of the phonons with wave numbers close to ±2kp 
only. Notin g th at the dis placement of ions Ui(x) is given 
from Eqs. (O) and (O) by, 



Ui(x) 
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{Ai(x)t 



A*(x)e"^}, (2.10) 



we obtain the ionic contribution to the charge density 
(2.4) 

modulation as, 



pS? = p, 



PQi - -P 



dui 



dx 
P^qQ 

g 



{Ai(x)e iQx + A*(x)e- iQx } 



(2.11) 



(2.12) 



where p is the average ionic charge densit y. No te that we 
neglected the derivative d x Ai (x) in Eq. ( 2.12 ), since its 
effect is negligible as compared to QAi. 

As a whole, the total charge density poi and Pqi are 
given by the summation of the electronic and the ionic 
contributions as, 



Poi{x) =p-e {R* <J R ia + L* a Li. 
p Qi = A*(x)e- iQx +Ai(x)e u 



(2.13) 
(2.14) 



where 



B. Interchain Coupling 

1. Coulomb Interaction 

As a source of the interchain coupling we first consider 
the Coulomb interaction @,|H. The Coulomb interac- 
tion is introduced by adding the following action, 

+ |dr|dri-|V^(r,r)| 2 , (2.8) 



Ai(x) = r)Ai(x) - eL* a R ia , 



(2.15) 



with 77 being pctQQ/g. 

In order to clarify the nature of the Coulomb interac- 
tion, first we trace out <f(x, r). One can easily see that 
the following term is generated, 



S c = - / drJ2 dx dx 



y Jo Jo 



j Pi(x)Pj(x') 
(x - x') 2 + dfj 

(2.16) 



where dy is the distance between the i-th and the j-th 
chain and the summation over i and j should be taken 
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over all the chains. This expression can be rewritten 
using Fourier transformation as, 

dr 



H l H ^3 (q)pi (<?)pj ■(-<?) » 



(2.17) 



Pi{q) = / dxpi(x) e lqx , 



;(?) = 



dx 
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= 2if (d y M) (L— foo), 



(2.18) 



(2.19) 



where -Kq(x) is the modified Bessel function. We assume 
that da = do ^ in or der to avoid divergence which oc- 
curs when i = j in Eq. ( 2. IE ) . The divergence arises from 
the fact that we considered the purely one-dimensional 
case, i.e., infinitely thin chains. In reality, do corresponds 
to the actual thickness of the chains, which may be given 
by the size of ions (or, more precisely, the extension of 
the electronic wave function on an ion). 

By subst itutin g pi(x) = poi(x) + PQi{x) into Eqs. 



( |2.17|) and fl2.18|) , we obtain 
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dr 



^ L ^ 



-2v ij (Q + q)A*(q)A j (q) 
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where poi{q) and Ai(q) are the Fourier components of 
Poi(x) and Ai(x), respectively. We have neglected the 
cross term between poi (x) and PQi (x) , since Ri a and Li a 
are slowly varying. 



Note that the interactions in S, 



(i) 



and Sq^ have differ- 



ent range in dij . As is shown in Eq. ( 2.19|) , Vij (q) is given 
by 2K (dij\q\) which decays exponentially for dij\q\ > 1. 

Since the wave number q in Sq can take the value q ~ 0, 

( 2) 

^o(dy|g|) is not negligible even for large dij. In Sq , 
on the other hand, Vij(Q + q) equals approximately to 
2_ftT (Qdij ) and this is negligible for d^ > 1/Q = 1/ (2kp). 
The Fermi wavelength is usually the order of the lattice 
constant and we can limit the summation with respect 
to i and j in 5q to adjacent chains only. 

We rewrite the long range part Sq\ restoring the 
scalar potential, as 



/' " pi 

/ dT ^ 
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dx 



■i<Pi {p-e (R* a R t a + L*a L ia)} 



{ipi +a - <Pif 

dl 



(2.21) 



where a = {y, z} and i + y and i + z symbolically denote 
the chains shifted by one lattice vector in positive y- and 



z-direction from the i-th chain, respectively. Note that 
only the slowly varying part of the scalar potential with 
a wave vector much smaller than kp should be taken into 
account here. 

(2) 

The short range part Sq can be divided into two 
kinds of terms: intrachain term (j — i) and the nearest- 
neighbor term (j = i + a). The intrachain term has three 
parts: ion-ion interaction, electron-ion interaction and 
electron-electron interaction. The first and the second 
one are absorbed into the phonon energy Tiujq and the 
electron-phonon coupling constant g, respectively. In this 
paper we assume that the electron-electron interaction is 
small and negligible. Hence we suppress the intrachain 
term in the following. 

The nearest-neighbor term (j = i + a) can be written, 
using the approximation, 
as, 



, (Q + q) — > v ir+a (Q)=v 



4 2) = J o dr Y, d <* j dx { A < ( x ) A *+<* i x ) 



J 4| +Q (x)A(x)}. 



2. Electron hopping 



(2.22) 



Next we discuss the effect of electron hopping between 
chains. We consider only the hopping between adjacent 
chains and introduce the following term, 



r 0h pL 
Shop = j dr / dx 



{t a (R* +a<7 Ri„ + L* +aa L i<T ) + c.c.} , (2.23) 

where t a is transfer integral between the i-th and the (i + 
a)-th chains and c.c. stands for the complex conjugate. 



C. Chiral Transformation 

In integrating out the electronic degree of freedom, it 
is convenient to employ chiral transformation so that the 
phase 0i(x) is eliminated from Aj(x). In this paper we 
study the temperature region not too close to T c . Hence 
the amplitude of the order parameter can be treated as 
a constant. We write Aj(x) = Aj e 19 ^^ in the following. 
The chiral transformation is introduced as follows: 

R ia — ► exp(i0 i /2)i? i(T , L ia — ► exp(-i6» i /2)L i(7 . (2.24) 

By this transformation S e - P is changed to, 
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s e . p = dr y2 dx 
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where 

Ki(x,T) 



*&° + hi 



km ^ + hi 



= hd T =F ihv F d x , 



h i 1 

hi = 7— i fep^^ + - (a^) 2 

2m 4 



(2.26) 

(2.27) 
(2.28) 



and vp = hkpjm is the Fermi velocity. 

In Eq. (2.27) we have neglected the term proportional 
to c* 2 .. As pointed out by Ishikawa and Takayama fjlf l 
this term is necessary to obtain the free energy which de- 
scribes electromagnetic response of the system correctly. 
On the other hand, if one neglects this term and ap- 
proximates the dispersion by a linear one, so-called chi- 

1<||. In 
lowing 



ral anomaly term must be taken into account jl^ 
this paper, we take the latter way and add the fo 
term to our action, 



S c h= dT J2 dx —d m 9i(x)<pi{x). (2.29) 



The terms Sq and Shop are transformed as, 



S, 



(2) 



0H pL 

At > v„ I dx 

2rj 2 AiA l+a cos(9 l+a - 9i) 

~ r)e{Ai {R* +aa L l+acr e- t{ 

+ (i < — ► i + a) I 



-0.) 



C.C. 



+ C.C 



} 



(2.30) 



,/3ft ,L 

Shop 



e i+a -9i 



1, a 
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+e* 2 L* i+aa Lia) +c.c.|. 



(2.31) 



III. DERIVATION OF THE GL FREE ENERGY 

The effective free energy of the condensate can be ob- 
tained by integrating out the electronic degree of freedom 
in the following way, 

F[0 u <pi, Ai] =-iln J Y[v[^l^ ia }e- s ^/\ (3.1) 

where Stot is the total action. In performing this path in- 
tegral, the following perturbative treatment is available, 



Z{fi) = j V[*l, exp [-{S + SS}/h] 



n=0 



(-irss\ _ Sa/R 



i\h n 



MP) x J2 



(-l)"(*s") t 

n\h n 



OO 

Z ((3)xexp[-J2(SS r - 



Z (p) xexp[W{^,^,AJ] , 



(3.2) 



where Sq and SS are the unperturbed and perturbed 
part, respectively. The bracket (••■)„ indicates that only 
the "connected diagrams" should be taken into account 
and Zq(/3) is the unperturbed partition function. From 
this formula, the effective free energy F is given by 
F = F Q — W/(3 where F = - ln(Z (/?))//?. From now on 
we assume that the amplitude of the order parameter Aj 
does not depend on the chain index and write Aj = A. 

Based on the above mentioned method, the effective 
free energy can be derived. Here, for simplicity, we dis- 
cuss only the resulting free energy, leaving the details of 
the calculations for Appendix [a|. The free energy consists 
of two parts, the intrachain part -Fintra and the interchain 
part -Pinter- The intrachain part is given by, 



F 



it 1 1 fa — ^ 



dx 



hvp 2 



e 2 (l - / s ) 2 ief s 

H r & H o x 6npi 

irnvp ir 

where the condensate fraction / s behaves as, 



fs 



1 

in 2 
1 - 



/3 2 A 2 C(3,i) 
V / 27r/3Aexp(- 



■0A) 



for PA < 
for PA > 



(3.3) 



(3.4) 



with C(3, 1/2) being the zeta function. (See Appendix [B] 
for details.) The interchain part -Pinter is given by, 

-Pinter = const. + ^ J dx KT a cos(0 i+a - 9i), (3.5) 



where r Q is the interchain coupling constant given by, 



r„ 



2v a A 2 ( poiqQ ehujc 



V 9 



2ff 2 



j-2 J8- 

irn vp 



(3.6) 



As one can see from Eqs. (3.5) and (3.6), the CDW 
prefers the phase difference of the neighboring chains to 
be 7r since r a is positive. In this paper we assume that 
the chains form a bipartite crystal and shift the phase of 
the chains on one sublattice by n. In this new notation 
T Q is changed to — r a , leading to the following expression 

for .Pinter, 



4 



-Pinter = const. — / dx flT a 

i,a 



cos 



'i-\-OL 



(3.7) 



where F^ eI is the change of the condensation energy due 
to interchain coupling. 

Next we take the continuum limit with respect to the 
chain index i by writing 6(r) instead of 9i{x). The total 
free energy becomes, 



dr 



K 



+ieJpd x e + — {|V^| 2 + (A - 2 + A^V 2 } 



(3.8) 



where K = hv F N ± f s /(2ir) and J = N±f s /-n with N ± 
being the cross sectional density of the chains. The 
anisotropy parameters j y and j z are given by K 7 2 = 
hTad^ (a = y,z). Here we have included the energy 
of electric field, i.e., the last two terms of Eq. (2.21) 
and also introduced the external charge p ext (r) which is 
a source of the external electric field. 

The first and the second term of T a in Eq. (3J3) are 
proportional to A 2 and / s , respectively. As is known for 
the case of superconductivity, these quantities, namely 
the square of the gap and the superfluid density, have dif- 
ferent temperature dependence at T <C T c , although the 
relation / s oc A 2 holds near T c |l7|]. Therefore 7 2 , which 
is proportional to T a / f s , has a weak temperature depen- 
dence at T <C T c arising from the difference between A 2 
and / s . This temperature dependence, however, is al- 
most negligible because the temperature dependence of 
A 2 and / s are both exponentially suppressed at low tem- 
peratures. 

In Eq. (3.8) Ao is the electrostatic screening length by 
the quasiparticles, which is given by Aq 2 = A^p(l — / s ) 
where A^p = 8e 2 N±/(hvF ) is the Thomas-Fermi screen- 
ing length of the individual 1-D chain. At T = 0, / s 
tends to 1 and Ao diverges. This is because we have as- 
sumed the complete nesting in our calculation. In actual 
CDW's, there exists residual carriers which does not con- 
tribute to CDW ordering even at OK. Such materials are 
called "metallic" CDW's. In this case, the divergence of 
Ao is suppressed. We have introduced A rcs to incorporate 
the contribution of these carriers. 

The amplitude fluctuation of the order parameter can 
be treated in a similar way as we have treated the phase 
fluctuation in this paper. For purely one-dimensional sys- 
tems this was carried out by Brazovskii and Dzyaloshin- 
skii |l8j . The extension of the present calculation to the 
case with the amplitude mode may be straightforward. 



IV. SUMMARY 

In this paper we have derived the effective free en- 
ergy of the phase of CDW employing the path integral 
method. Especially, a careful treatment of the interchain 
coupling, i.e., Coulomb interaction and electron hopping, 
was presented. We have clarified the followings: 

1. The Coulomb interaction is divided into two parts, 
which correspond to slow modulation (q ~ 0) and 
fast oscillation (q ~ 2kp) of charge density. The 
former determines the large scale electromagnetic 
response of the CDW's and the latter plays an im- 
portant role in the elastic coupling between CDW's 
on adjacent chains. 

2. The electron hopping between neighboring chains 
also contributes to the elastic interchain coupling 
of CDW. 

3. The contribution of the Coulomb interaction and 
the electron hopping to the interchain coupling 
have slightly different temperature dependence at 
low temperatures. As a result the anisotropy pa- 
rameters also become weakly temperature depen- 
dent at low temperatures, although its dependence 
is negligibly small. 

4. The quasiparticles give a finite electrostatic screen- 
ing length. The length is given by the Thomas- 
Fermi length above T c but is enhanced in lower tem- 
peratures due to the suppression of the density of 
states by the energy gap of CDW. In "semiconduct- 
ing" CDW's, the quasiparticle screening length is 
divergent at T = but it remains finite in "metal- 
lic" ones. 
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APPENDIX A: INTEGRATION OVER 
ELECTRONIC FIELD 



From Eqs. (|2~2l| ), ( ^25|) , (|2~29| ), ( ^30|) and (|2~3l| ), we 

obtain the total action of the system. For simplicity we 
set Ai — ► A in the following. We divide the total action 
into two parts, unperturbed part Sq and perturbe d par t 
SS. Here we assign Sq the part of Se-p in Eq. ( 2.25| ) 
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which does not include Oi. The rest of the terms are 
assigned to SS. 

The unperturbed action is given as follows, 



So f I 
— = drdr' > / dxdx 1 

& Jo V Jo 



*U*> r) ■ {~G(x -x'.r- t')}- 1 ■ Vurtfy), 



2 fc,e n 



where, 



G(k,ie n ) 



G R (k,ie n ) F(k,ie n ) 
fi(k,ie n ) G L {k,ie n ) 



and the Green functions are given by 
-ft (ihe n + f) 



G R (k,ie n ) 
G L (k,ie n ) 
F{k,ie n ) 



-ft (ihe n - £) 
(fe„) 2 + ^(fc) 2 ^ 

-TiA 
(fe„) 2 + £(fc) 2 ' 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



with J5(fc) = V^ + A 2 , £ = ^w F fc and s n = (2n 

l)7T/(^). 

The perturbed term 5S 1 is given as follows, 



= Si + 5*2 + S, 



(2) 



Si 



lop, 



51 = / dT E / dzXiW (i^iiia + ^ii.) , (A6) 

i 

5 2 = -i[dr^[dx & [p - ^d x 6i] , (A7) 



where Xi( x ) i s defined by 



Xi{x) = —ipi 



-vf9 x 6i + - — (d x 8i) 2 . 
2 8m y ' 



(A8) 



S ( c } and Shop are given by Eq. ( |2.30| ) and Eq. (^3lj) 



respectively. The energ y of the electric field, i.e., the last 
two terms of Eq. (2.21) is not included here. 

After the integration with respect to the electronic de- 
gree of freedom we obtain the following results: 



(5i) c = J dT Yl [ dxn eXi(x), 



(A9) 



where n c is the equilibrium electron density on a chain 
which is assumed to be a constant independent of i and 
given by 2kp/ix. In Eq. (|A9|), the first term of Xi( x ) 
should be cancelled by the ionic background charge, i.e., 
the first term of Eq. (A7) for charge neutrality. In the 
present calculations, the cancellation of second term of 



Xi(x), which is proportional to d x &i, is not complete. 
This term should vanish in order to maintain the sta- 
bility of CDW with a wave number 2kp. We consider 
that this cancellation is fulfilled if we take into account 



the <9 2 -term which we have neglected in Eq. (2.26). 

The second order term with respect to Si is given as 
follows, 

({SiH c = E/ dXdX' Xl (X) Xt (X') 

i 

x2^-g R (x - x')g R (x' - x) 

-G L {X-X')G L {X' -X) 

-2T{X - X')T{X' - X)\, (A10) 

where we have introduced X = (x, t) for simplicity. Since 
Xi( x ) is slowly varying as compared to G R , G L or T , we 
can use the following treatment, 

g R (x ~ x')g R (x' - x) 

-^S(X-X')x f dX G R (X)G R (-X) 



It i, 



(All) 



This treatment corresponds to GL expansion. As we see 
in the next section, the following quantities are rewritten 
in terms of the the condensate fraction / s as, 



(SftL 



L.J2{s R (k,is n ) 2 +G L {k,ie n f} 



k,Sn 



1 

ttvf 



1 - 



fifth 



(A12) 
(A13) 



which yields the following expression, 

= /dr^|dx|-/ B ^(W 2 
-f s —p<Pid x 9i - (1 - / s ) -^ipj 

KIT, TTVFft 

= Wi. 



(A14) 



The terms which includes two chains are given as fol- 
lows, 



,(2) 



fin 



dr / dx 



f pcxqQ eftiOQ 



\ 9 



25 2 



x2A z cosf 



(A15) 
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where we have used the relation, Eq. (2.7). We also 
obtain for the hopping term, 



({W 2 ) c =4/ dXcLY']rM 2 { 



~g R {x - x')g R {x' - x) 

—Q L {X — X')Q L (X' — X) 

-T*{X - X')T{X' - X)e~ l( - e '+"- 6 ^ 

-Fi(x - x')t*(x' - xy {8 ^- 6 * ] } 

*/«Ew 2 {^(i-f) 



2t:vi 



■ cos 



,)}. 



(A16) 



Collecting these terms we obtain the interchain terms 



as, 



2|t„| 2 ( l _f 1 



h z irvF V ^ 

r Q COs(9 i+a 



Wo 



where 



2v a A 2 f pa Q Q eTiojQ 



K\ 2 



(A17) 



-U (A18) 



h V g 2g 2 J Trh 2 v F ' 

Here we assumed Oi and tpi to be classical variables and 
neglected their r-dependence. 

The intrachain part -Fi n t ia and the interchain part -Pinter 
of the free energy are given, respectively, by i<i n t r a = 
-Wi/P and F intor = -W 2 /0. 



APPENDIX B: CALCULATION OF THE GL 
COEFFICIENTS 



In this section we derive Eq. (|A12|) and (|A1^) . We 
introduce a cut off energy, hvp\k\ ~ A, which should be 
set to in the end of the calculation. 

Eq. ( A12| ) can be obtained in the following way, 



1 —J2{G R (k,is n ) 2 +g L (k,ie n ) 2 } 

-(hs n ) 2 +e 



(3HL 
2 



E— / d ^ 

■~ KVF Jo 



= h + h 



{(hs n y + E(k) 2 } 2 ' 
■ e (hen? - A 2 

_{{he n f + E(k) 2 } 2 
A 2 

{(he n ) 2 + E(k) 2 } 2 



(Bl) 



-A 



I\ is calculated as, 

1 ' pli^TTVF Va 2 + {he n ) 2 + A 2 
-1 A 



: tanh 



/VA 2 + A 2 



irv F VA 2 + A 2 
— (A — >oo), 

TTVf 



using the identity, 

00 

tanh z — z 



1 



{(2n- 1)tt/2} 2 + z 2 



(B2) 



(B3) 



I2 is also calculated in a similar way as, 
A 2 ^ 2h 

1 ~Jm^^F~ 



2{(he n ) 2 + A 2 }{(he n ) 2 + A 2 + A 2 } 

1 A 

arctan 



2^/(he n ) 2 +A 2 

A 2 1/1 J3A 

r t tanh 

2tr; f AVA 2 

1 



tanh 



v/(fe„) 2 + A 2 



/VA 2 + A 2 



+ 



A 2 

itvf/3 



VA 2 + A 2 2 
£{(fe„) 2 + A 2 }-i 

A 



x arctan ■ 



1 A 2 

TTVp (3 

x arctan 



V(^n) 2 +A 2 ' 

^{(fe„) 2 + A 2 }-l 

En 

A 



V(fe„) 2 + A 2 
We define the condensate fraction / s by 

/ s = 2irv F h 

= 4 V f\ ]r A 2 

P^Jo Ue + (h£n) 2 +A 2 } 2: 



(A — ► 00). (B4) 



^£{(fe„) 2 + A 2 }-i 



x arctan ■ 



v/(fe„) 2 + A 2 



(B5) 



Then from Eq. flB2|) and we obtain Eq. ( |A12[ ) 
In a similar way we obtain, 



0hL 



(B6) 
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which reduces to Eq. (A13). 

Next we examine the behavior of / s in two limiting 
cases, /3A < 1 and /3A > 1. 



1. The case of /J A < 1 



In this case, we can simply expand Eq. (B5) with 
respect to A and obtain, 



ttA 2 



/3 ^ (he n r 

4tT 2 U S j 



Q9A) 

4tt 2 



2 00 



ra=l 1 2^ 



(B7) 



2. The case of /?A > 1 

In this case the following method is more useful |l9| 
From Eq. (lB5|), we obtain, 



A 2 



cA 



U 2 + (fen) 2 +A 2 } 2 ' 

x 9 /- A , 1 , /IE 



/3 2 (£ 2 + A 2 ) 



-A- 



_9_ 
"<9A 



^ A d^|l + 2E(-l) m e^"j 



log 



A + VA 2 + A 2 



+2 £(-!)' 



Va 2 +a 2 /a 



de- 



-/3Aem 



OO 

l + 2/?A^(-l) m mKi(/3Am) (A — ► oo) 
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27r/3Ae 



-f3A 



(B8) 



where E = v^ 2 + A 2 and e = E / A. K x (z) is the modi- 
fied Bessel function given by, 



dz 
d 



K {z) 



-zy 



dzj 1 ^y 2 -l 



(B9) 



In the last line of Eq. ( p38| ) , we have employed the asymp- 
totic form of Kx(z). 
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